Learning apparatus and method for learning a model corresponding to real number time-series input data

ABSTRACT

An apparatus, a computer readable medium, and a learning method for learning a model corresponding to time-series input data, including acquiring the time-series input data, which is a time series of input data including a plurality of input values, propagating, to a plurality of nodes in a model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point, calculating a node value of a first node among the plurality of nodes by using each propagated value propagated to the first node, and updating a weight parameter used to calculate each propagation value propagated to the first node, by using a corresponding input value and a calculated error of the node value at the one time point.

BACKGROUND 1. Technical Field

The present invention relates to an apparatus and a method for learning a model corresponding to real number time-series input data.

2. Related Art

A neural network, such as a Boltzmann machine and the like, can learn a model corresponding to data input in time series. In particular, a dynamic Boltzmann machine is expected to be able to achieve high learning power through machine learning.

Such a Boltzmann machine or the like handles binary data with values of 1 or 0 based on probability, for example, and therefore it is difficult for such a Boltzmann machine or the like to be adapted to real number time-series data.

SUMMARY

According to an embodiment of the present invention, a learning method for learning a model corresponding to time-series input data is provided, the method comprising acquiring the time-series input data, which is a time series of input data including a plurality of input values, propagating, to a plurality of nodes in a model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point, calculating a node value of a first node among the plurality of nodes by using each propagated value propagated to the first node, and updating a weight parameter used to calculate each propagation value propagated to the first node, by using a corresponding input value and a calculated error of the node value at the one time point. Also provided are an apparatus and a computer readable storage medium.

According to another embodiment of the present invention, a learning method using a learned model corresponding to time-series input data of a learning target is provided, the method comprising acquiring time-series input data of an investigation target, which is a time series of input data including a plurality of input values, propagating, to a plurality of nodes in the model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point; calculating a node value of a first node among the plurality of nodes by using each propagated value propagated to the first node, and detecting peculiarity of an input value corresponding to the first node and input at the one time point, by using the node value of the first node and a variance parameter for indicating a variance in a probability distribution of input values learned using the time-series data of the learning target. Also provided are an apparatus and a computer readable storage medium.

The summary clause does not necessarily describe all features of the embodiments of the present invention. The present invention may also be a sub-combination of the features described above.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows an exemplary configuration of a learning apparatus according to an embodiment;

FIG. 2 shows an exemplary configuration of a model according to an embodiment;

FIG. 3 shows a flow of an operation of the learning apparatus according to an embodiment;

FIG. 4 shows an example of structures of time-series data and training data for use in learning in an embodiment;

FIG. 5 shows a first modification of the learning apparatus according to an embodiment;

FIG. 6 shows a modification of the model according to an embodiment;

FIG. 7 shows an example of a temporal change in the update parameter γ_(i,l) according to an embodiment;

FIG. 8 shows an example of a temporal change in the update parameter α_(i,j,k), according to an embodiment;

FIG. 9 shows an operational flow of a learning apparatus of a modification according to an embodiment; and

FIG. 10 shows an exemplary hardware configuration of a computer according to an embodiment of the present invention.

DESCRIPTION OF EXEMPLARY EMBODIMENTS

Hereinafter, some embodiments of the present invention will be described. The embodiments do not limit the invention according to the claims, and all the combinations of the features described in the embodiments are not necessarily essential to means provided by aspects of the invention.

FIG. 1 shows an exemplary configuration of a learning apparatus 100 according to an embodiment. The learning apparatus 100 may be an apparatus for learning a model 10 corresponding to time-series input data. The learning apparatus 100 may be configured to learn a model based on a Boltzmann machine by supplying time-series data to nodes of the model. The learning apparatus 100 includes an acquiring section 110, a supplying section 120, a storage section 130, a computing section 140, and a learning processing section 150.

The acquiring section 110 may be configured to acquire time-series input data. Time-series input data may be, for example, a data sequence in which a plurality of pieces of data are arranged along a time axis, such as moving image data. The acquiring section 110 may be connected to a device operated by a user or a device (e.g., sensor) that detects and/or outputs time-series data and may acquire the time-series input data from such a device. Alternatively, the acquiring section 110 may read and acquire time-series input data stored in a storage device in a predetermined format. Alternatively, the acquiring section 110 may be connected to a network and acquire time-series input data via the network. The acquiring section 110 may also store the acquired time-series input data in a storage device included in the learning apparatus 100.

The supplying section 120 may be configured to supply a plurality of input values corresponding to input data at one time point in the time-series input data to a plurality of nodes of a model. The supplying section 120 is connected to the acquiring section 110 and may handle, as training data, input data at one time point in the received time-series input data and supply input values at the one time point to corresponding nodes of the model. Input data at one time point may be the temporally newest data in a training data set for use in learning. Alternatively, input data at one time point may be temporally intermediate data in a training data set for use in learning. That is, input data at one time point may be selected arbitrarily from the time-series data.

The storage section 130 may be configured to store values of hidden nodes of the model in correspondence with a plurality of time points in the time-series input data. The storage section 130 may sample the values of hidden nodes corresponding to one time point and store these values in the hidden nodes corresponding to this time point. The storage section 130 may store the sampled values respectively in the hidden nodes corresponding to the time points for each time point.

The computing section 140 may be configured to compute a conditional probability of each input value at one time point on a condition that an input data sequence has occurred. Here, in the time-series input data, the pieces of data input to the respective nodes of the model at one time point are referred to as input values at one time point, and the pieces of data input to the respective nodes at the time points before the one time point are referred to as the input data sequence. A model used by the learning apparatus 100 may have a weight parameter between (i) a plurality of hidden nodes and a plurality of input values corresponding to input data at each time point prior to the one time point in an input data sequence and (ii) a plurality of hidden nodes corresponding to the one time point and a plurality of input nodes.

The computing section 140 may be configured to compute a conditional probability of each input value at one time point, on the basis of an input data sequence before the one time point in the time-series input data, the stored values of hidden nodes, and the weight parameter of the model. Furthermore, the computing section 140 may be configured to compute a conditional probability of the value of each hidden node at one time point on a condition that an input data sequences has occurred, based on an input data sequence before the one time point in the time-series input data and the weight parameter of the model.

The learning processing section 150 may be configured to increase a conditional probability of input data at one time point occurring on a condition that the input data sequence has occurred, by adjusting the weight parameter of the model. The learning processing section 150 may further adjust bias parameters which are given respectively to the plurality of nodes and hidden nodes of the model. The learning processing section 150 may supply the adjusted weight parameter and bias parameters of the model to a storage device, such as an external database 1000, to store these parameters in the storage device.

The above-described learning apparatus 100 according to the present embodiment may be configured to learn the model by adjusting the weight parameter and bias parameters of the model, based on input data at one time point in the time-series input data. The model, according to the present embodiment, is described with reference to FIG. 2.

FIG. 2 shows an exemplary configuration of a model 10 according to the present embodiment. The model 10 includes a plurality of common layers 12. FIG. 2 shows an example including a total of T common layers 12. The model 10 may include a finite number of common layers 12. Each common layer 12 includes an input layer 14 and a hidden layer 16.

Each input layer 14 may be a layer corresponding to the time-series data. Each input layer 14 may correspond to a respective time point in the time-series data. Each input layer 14 may include a predetermined number of nodes. For example, the 0-th input layer may be a layer corresponding to input data at one time point in the time-series data. The 0-th input layer may include a plurality of nodes corresponding to the number of input values in this input data.

A total of T−1 input layers 14 other than the 0-th input layer 14 among the plurality of input layers may be input layers 14 corresponding to the input data sequence before the one time point in the time-series input data. For example, the −1st input layer may correspond to input data at a time point that temporally precedes the one time point by one time point, and the (−δ)-th input layer may correspond to input data at a time point that temporally precedes the one time point by δ time points. That is, a total of T−1 input layers other than the 0-th input layer each have the same number of nodes as the 0-th input layer and are respectively supplied with input values of corresponding input data values in the input data sequence, for example.

Each hidden layer 16 may correspond to a respective time point in the time-series data. For example, the 0-th hidden layer may be a layer corresponding to one time point in the time-series data. FIG. 2 shows an example including a total of T hidden layers 16. Each hidden layer 16 may include one or more hidden nodes, and the storage section 130 may store the values sampled at the one time point.

A total of T−1 hidden layers other than the 0-th hidden layer among the plurality of hidden layers 16 may be hidden layers 16 corresponding to time points before the one time point in the time-series data. For example, the −1st hidden layer corresponds to a time point that temporally precedes the input data of the one time point by one time point, and the storage section 130 stores the values sampled at the time point that temporally precedes the one time point by one time point. Furthermore, the (−δ)-th hidden layer may correspond to a time point that temporally precedes the input data of the one time point by δ time points, and the storage section 130 may store the values sampled at the time point that temporally precedes the one time point by δ time points. That is, a total of T−1 hidden layers other than the 0-th hidden layer each have the same number of nodes as the 0-th hidden layer and are respectively supplied with values of corresponding hidden nodes, for example.

As an example, in the case where the time-series input data is moving image data, the last image data of the moving image data corresponds to the 0-th input layer, and a plurality of nodes of the 0-th input layer each receive corresponding pixel data of the image data. Furthermore, the 0-th hidden layer corresponds to the final time point of the moving image data, and the storage section 130 may store values sampled at this final time point in the hidden nodes of the 0-th hidden layer.

In addition, the −1st input layer is supplied with image data that immediately precedes the last image data, and a plurality of nodes of the −1st input layer each receive corresponding pixel data of the immediately preceding image data. Furthermore, the −1st hidden layer corresponds to the time point that immediately precedes the final time point, and for each of the plurality of nodes of the −1st hidden layer, the storage section 130 may store the values sampled at this immediately preceding time point. Similarly, the plurality of nodes of the (−δ)-th input layer each receive corresponding pixel data of image data that precedes the last image data by δ images, and the plurality of nodes of the (−δ)-th hidden layer each store corresponding sampling values at the time point that precedes the last time point by δ time points.

FIG. 2 shows an example in which each common layer 12 includes an input layer 14 and a hidden layer 16, but instead, one or more common layers 12 need not include a hidden layer 16. In such a case, the 0-th common layer to the (−m)-th common layer includes input layers 14 and hidden layers 16, and the (−m−1)-th common layer to (−T+1)-th common layer may include input layers 14.

The plurality of nodes in the 0-th input layer 14 and/or the plurality of hidden nodes in the 0-th hidden layer 16 may each have a bias parameter. For example, the j-th node j in the common layer 12 has a bias parameter b_(j).

The plurality of nodes in the 0-th input layer 14 and the nodes of the hidden layer 16 corresponding to the input data sequence and layers corresponding to the input data sequence before the one time point may respectively have weight parameters therebetween. There need not be weight parameters between the plurality of nodes in each input layer 14 and hidden layer 16.

Similarly, the plurality of nodes in the 0-th hidden layer 16 and the nodes of the hidden layer 16 corresponding to the input data sequence and layers corresponding to the input data sequence before the one time point may respectively have weight parameters therebetween. That is, the plurality of nodes of the 0-th common layer 12 and the nodes of the plurality of common layers 12 before the one time point may respectively have weight parameters therebetween.

FIG. 2 shows a concept of a weight parameter W_(ij) ^([δ]) between the node j of the 0-th input layer and a node i of the (−δ)-th layer. FIG. 2 shows an example in which the model 10 has the same number of input layers 14 and layers 16, each input layer 14 includes I nodes, and each hidden layer 16 includes H hidden nodes. In the present embodiment, the input layers 14 and hidden layers 16 are expressed by one common layer 12 that has a plurality of nodes x_(j) ^([t]). The first to I-th nodes (1≤j≤I) of the common layer 12 indicate the nodes of the input layer 14, and the (I+1)-th to (I+H)-th nodes (I+1, j, I+H) indicate hidden nodes.

For example, the weight parameter W_(ij) ^([δ]) shown in FIG. 2 indicates a weight parameter between two nodes of an input layer 14. Similarly, the weight parameter W_(ij) ^([δ]) may include weight parameters from an input node to a hidden node, from a hidden node to an input node, and between to hidden nodes. A specific example of a weight parameter W_(ij) ^([δ]) is shown in Expression 1 below. In the present embodiment, a weight parameter from an input node or hidden node to an input node is referred to as a “weight parameter to an input node,” and a weight parameter from an input node or hidden node to a hidden node is referred to as a ‘weight parameter to a hidden node.”

$\begin{matrix} {{W_{ij}^{\lbrack\delta\rbrack} = {{\hat{W}}_{ij}^{\lbrack\delta\rbrack} + {\hat{W}}_{ij}^{\lbrack{- \delta}\rbrack}}}{{\hat{W}}_{ij}^{\lbrack\delta\rbrack} = \left\{ \begin{matrix} 0 & {{{if}\mspace{14mu}\delta} = 0} \\ {\sum\limits_{k \in K}\;{u_{i,j,k}\lambda_{k}^{\delta - d_{ij}}}} & {{{if}\mspace{14mu}\delta} \geq d_{ij}} \\ {\sum\limits_{l \in L}\;{{- v_{i,j,l}}\mu_{k}^{- \delta}}} & {otherwise} \end{matrix} \right.}} & {{Expression}\mspace{11mu} 1} \end{matrix}$

Here, u_(i,j,k) and v_(i,j,l) are learning parameters that are learning targets, for example. Furthermore, λ_(k) ^(t1) and μ_(l) ^(t2) are predefined parameters that change in a predetermined manner in accordance with a time point difference δ between the hidden nodes and input data in the input data sequence before the one time point and the hidden nodes and input data at the one time point (t1=δ−d_(ij), t2=−δ). That is, the weight parameter W_(ij) ^([δ]) may be a parameter based on the learning parameters u_(i,j,k) and v_(i,j,l) and the predefined parameters λ_(k) ^(t1) and μ_(l) ^(t2).

The weight parameter W_(ij) ^([δ]) may be a parameter based on a positive value, which is based on a product of the first learning parameter u_(i,j,k) and the first predefined parameter λ_(k) ^(t1), and a negative value, which is based on a product of the second learning parameter v_(i,j,l) and a second predefined parameter μ_(l) ^(t2). Specifically, in the case where the time point difference δ is greater than or equal to a predetermined delay constant d_(ij), the weight parameter W_(ij) ^([δ]) may be a positive value based on a product of the first learning parameter u_(i,j,k) and the first predefined parameter λ_(k) ^(t1). In the case where the time point difference δ is less than the delay constant d_(ij) and is not equal to 0, the weight parameter W_(ij) ^([δ]) may be a negative value based on a product of the second learning parameter v_(i,j,l) and the second predefined parameter μ_(l) ^(t2). In addition, in the case where the time point difference δ is equal to 0, the weight parameter W_(ij) ^([δ]) may be equal to 0.

In addition, in the case where the time point difference δ is greater than or equal to the predetermined delay constant d_(ij), the weight parameter W_(ij) ^([δ]) may be based on a plurality of positive values that are based on the products u_(i,j,k)·λ_(k) ^(t1) of a plurality of sets of the first learning parameter u_(i,j,k) and the first predefined parameter λ_(k) ^(t1) respectively from among the plurality of first learning parameters u_(i,j,k) and the plurality of first predefined parameters λ_(k) ^(t1). In addition, in the case where the time point difference δ is less than the predetermined delay constant d_(ij) and is not equal to 0, the weight parameter W_(ij) ^([δ]) may be based on a plurality of negative values that are based on products v_(i,j,l)·μ_(l) ^(t2) of a plurality of sets of the second learning parameter v_(i,j,l) and the second predefined parameter μ_(l) ^(t2) respectively from among the plurality of second learning parameters v_(i,j,l) and the plurality of second predefined parameters μ_(l) ^(t2).

A predefined parameter may be a parameter based on a value obtained by raising a predetermined constant to the power of a value based on the time point difference δ. The first predefined parameter λ_(k) ^(t1) is a parameter whose value gradually decreases as the time point difference δ increases, for example. In this case, the first predefined parameter λ_(k) ^(t1) may be a value obtained by raising a first constant λ_(k), which is greater than 0 and less than 1, to the power of a value obtained by subtracting the predetermined delay constant d_(ij) from the time point difference δ (δ−d_(ij)=t1). In addition, the second predefined parameter μ_(l) ^(t2) may be a parameter whose value gradually decreases as the time point difference δ increases, for example. In this case, the second predefined parameter μ_(l) ^(t2) may be a value obtained by raising a second constant 1.1), which is greater than 0 and less than 1, to the power of a negative value of the time point difference δ (−δ=t2).

The above-described model 10 according to the present embodiment may be configured to form a Boltzmann machine. That is, the model 10 may be a Boltzmann machine to which time-series data is applied. The model 10 may be a Boltzmann machine that includes hidden layers 16 into which are input values differing from the time-series data, in addition to the input layers 14 into which the time-series data is input. The learning apparatus 100 according to an embodiment learns the model 10 by adjusting the learning parameters u_(i,j,k) and v_(i,j,l) and the bias parameter b_(j) while sampling and storing the values of the hidden nodes, by using, as training data, input data at one time point that is supplied to the 0-th input layer of the model 10. A learning operation of the learning apparatus 100 is described with reference to FIG. 3.

FIG. 3 shows a flow of an operation of the learning apparatus 100 according to the present embodiment. In the present embodiment, the learning apparatus 100 may be configured to learn the model 10 corresponding to time-series input data and determine the learning parameters u_(i,j,k) and v_(i,j,k) and the bias parameter b_(j), by executing the processing steps of S310 to S360. In an embodiment, first, an example is described in which the determination of the weight parameters to the hidden nodes and the weight parameters to the input nodes is performed by the learning apparatus 100 using substantially the same operation.

First, the acquiring section 110 may acquire time-series data (S310). The acquiring section 110 may acquire time-series data of a duration equivalent to a total of T layers from the 0-th layer to the (−T+1)-th layer of the model 10. The acquiring section 110 acquires, for example, T pieces of image data in time-series that form the moving image data.

Then, the supplying section 120 may supply a plurality of input values corresponding to the input data of the time-series input data at one time point to the plurality of input nodes of the 0-th input layer of the model 10 (S320). Here, x_([1,I]) ^([0])(=x_(j) ^([0]), 1≤j≤I) denotes input data supplied to the 0-th input layer.

The supplying section 120 supplies, for example, I input values x_(j) ^([0]) corresponding to input data x_([1,I]) ^([0]) of the time-series input data at the most recent time point to the corresponding nodes j of the 0-th input layer (1≤j≤I). For example, the supplying section 120 supplies I pieces of pixel data included in the last piece of image data of T pieces of image data arranged in time series to form the moving image data to I nodes of the 0-th input layer. The supplying section 120 may supply a value of 1 or 0 as the pixel data to each node of the 0-th input layer. If the duration of the time-series input data is shorter than T, the supplying section 120 may supply the data to a number of layers from the 0-th input layer corresponding to the length of the time series, and may supply a value of 0, for example, to the nodes of the rest of the layers.

Then, the supplying section 120 may supply a plurality of input values corresponding to the input data sequence before the one time point to the plurality of nodes included in respective layers from the −1st input layer to the (−T+1)-th input layer of the model 10. Here, let x_(j) ^((−T, −1]) denote input data supplied to layers from the −1st input layer to the (−T+1)-th input layer (1≤j≤1). The term (−T, −1] indicates layers from the (−T+1)-th layer to the −1st layer. That is, the input data x_(j) ^((−T, −1]) in the time-series data denotes a history up to the input data x_(j) ^([0]), for example.

Next, the storage section 130 samples the values of a plurality of hidden nodes corresponding to the one time point, and respectively stores these values in the corresponding plurality of hidden nodes of the 0-th hidden layer (S330). The storage section 130 may arbitrarily or randomly input values of 1 or 0. The storage section 130 stores H sampled values in the corresponding hidden nodes j of the 0-th hidden layer, for example (I+1≤j≤I+H).

The storage section 130 may store the values of the hidden nodes before the one time point respectively in a plurality of nodes in each of the corresponding hidden layers from the −1st hidden layer to the (−T+1)-th hidden layer of the model 10. Here, let, let x_(j) ^((−T, −1]) denote the values of the hidden nodes stored in the layers from the −1st hidden layer to the (−T+1)-th hidden layer (I+1≤j≤I+H). That is, the values x_(j) ^((−T, −1]) input to the nodes of each common layer 12 before the one time point denote a history up to the input values x_(j) ^([0]) input to the nodes of the 0-th common layer, for example (1≤j≤I+H).

Then, the computing section 140 may compute conditional probabilities of each input value x_(j) ^([0]) (1≤j≤1) of an input node at the one time point, based on the input values x_(j) ^((−T, −1]) (1≤j≤I+H) of the plurality of nodes of the (−T+1)-th common layer to the −1st common layer and the weight parameter W_(ij) ^([δ]) (S340). The computing section 140 computes a probability <x_(j) ^([0])>_(θ) of the input value x_(j) ^([0]) (1≤j≤1) of the j-th node of the 0-th input layer being equal to 1 by substituting 1 for x_(j) ^([0]) in the following expression, based on the history x^((−T, −1]) (1≤j≤I+H) of the plurality of nodes of the common layer 12.

$\begin{matrix} {\left\langle X_{j}^{\lbrack 0\rbrack} \right\rangle_{\theta} = {{p_{\theta,j}\left( {x_{j}^{\lbrack 0\rbrack}❘x^{({{- T},{- 1}}\rbrack}} \right)} = \frac{\exp\left( {{- \tau^{- 1}}{E_{\theta,j}\left( {x_{j}^{\lbrack 0\rbrack}❘x^{({{- T},{- 1}}\rbrack}} \right)}} \right)}{1 + {\exp\left( {{- \tau^{- 1}}{E_{\theta,j}\left( {x_{j}^{\lbrack 0\rbrack}❘x^{({{- T},{- 1}}\rbrack}} \right)}} \right)}}}} & {{Expression}\mspace{14mu} 2} \end{matrix}$

In an embodiment, an example is described in which the input value x_(j) ^([0]) of each node is binary, e.g., 1 or 0, but the value of the input value x_(j) ^([0]) is not limited to these values. Furthermore, in the step for computing the conditional probabilities of each input value x_(j) ^([0]) (1≤j≤I) of the 0-th input layer, the computing section 140 may compute the conditional probabilities of the values x_(j) ^([0]) (I+1≤j≤I+H) of the 0-th hidden layer.

Expression 2 is derived as a Boltzmann machine from a known probability formula. For example, θ denotes a set of parameters to be computed, and the formula θ=(b_(j), u_(i,j,k), v_(i,j,l)) is established. In addition, T may be a parameter that is dependent on a known “system temperature” of the Boltzmann machine, and may be preset by a user or the like. Also, E_(θ,j)(x_(j) ^([0])|x^(−(−T,−1])) of Expression 2 is computed by using the following expression.

$\begin{matrix} {{E_{\theta,j}\left( {x_{j}^{\lbrack 0\rbrack}❘x^{({{- T},{- 1}}\rbrack}} \right)} = {{{- b_{j}}x_{j}^{\lbrack 0\rbrack}} - {\sum\limits_{t = {- T}}^{- 1}\;{\left( x^{\lbrack t\rbrack} \right)^{\top}\left( x^{\lbrack t\rbrack} \right)W_{:{,j}}^{\lbrack{- t}\rbrack}x_{j}^{\lbrack 0\rbrack}}}}} & {{Expression}\mspace{14mu} 3} \end{matrix}$

Here, “T′” denotes a transpose, “:” denotes 1 to n in a case where n (=I+H) denotes the number of nodes, and “:,j” indicates extraction of the j-th column. That is, the second term on the right side of Expression 3 is denoted by the following expression, for example.

$\begin{matrix} {{\sum\limits_{t = {- T}}^{- 1}\;{\left( x^{\lbrack t\rbrack} \right)^{\top}\left( x^{\lbrack t\rbrack} \right)W_{:{,j}}^{\lbrack{- t}\rbrack}x_{j}^{\lbrack 0\rbrack}}} = {{\left( {x_{1}^{\lbrack t\rbrack},x_{2}^{\lbrack t\rbrack},\ldots\mspace{14mu},x_{n}^{\lbrack t\rbrack}} \right)\begin{pmatrix} W_{1,j}^{\lbrack{- t}\rbrack} \\ W_{2,j}^{\lbrack{- t}\rbrack} \\ \vdots \\ W_{n,j}^{\lbrack{- t}\rbrack} \end{pmatrix}x_{j}^{\lbrack 0\rbrack}} = {\sum\limits_{i = 1}^{N}\;{\left( {{\sum\limits_{k \in K}{u_{i,j,k}\alpha_{i,j,k}}} - {\sum\limits_{l \in L}\;{v_{i,j,l}\beta_{i,j,l}}} - {\sum\limits_{l \in L}\;{v_{j,i,l}\gamma_{i,l}}}} \right)x_{j}^{\lbrack 0\rbrack}}}}} & {{Expression}\mspace{14mu} 4} \end{matrix}$

Here, α_(i,j,k), β_(i,j,l), and γ_(i,l) are denoted by the following expressions.

$\begin{matrix} {\alpha_{i,j,k} \equiv {\sum\limits_{t = {- T}}^{- d_{ij}}{\lambda_{k}^{{- t} - d_{ij}}x_{i}^{\lbrack t\rbrack}}}} & {{Expression}\mspace{14mu} 5} \\ {\beta_{i,j,l} \equiv {\sum\limits_{t = {{- d_{ij}} + 1}}^{- 1}{\mu_{l}^{t}x_{i}^{\lbrack t\rbrack}}}} & {{Expression}\mspace{14mu} 6} \\ {\gamma_{i,l} \equiv {\sum\limits_{t = {- T}}^{- 1}{\mu_{l}^{- t}x_{i}^{\lbrack t\rbrack}}}} & {{Expression}\mspace{14mu} 7} \end{matrix}$

Accordingly, P_(θ,j)(1|x_(j) ^((−T,−1])) obtained by substituting 1 for x_(j) ^([0]) of Expression 2 can be computed from Expression 3 by substituting 1 for x_(j) ^([0]) in expression 5. Note that predetermined initial values (for example, 0) may be substituted for the parameter set θ=(b_(j), u_(i,j,k), v_(i,j,l)). In this way, the computing section 140 can compute a conditional probability <x_(j) ^([0])>_(θ) of each input value x_(j) ^([0]) at the one time point which is denoted by Expression 2.

Then, the learning processing section 150 may adjust the parameter set θ=(b_(j), u_(i,j,k), v_(i,j,l)) (S350). When adjusting the bias parameter b_(j), the learning processing section 150 may determine a direction of the change in the bias parameter b_(j) by using the following expression.

$\begin{matrix} {{\frac{\partial}{\partial b_{j}}\log\;{P_{\theta}\left( x_{\lbrack{1,I}\rbrack}^{\lbrack 0\rbrack} \middle| x^{({{- T},{- 1}})} \right)}} = {\tau^{- 1}\left( {x_{j}^{\lbrack 0\rbrack} - \left\langle X_{j}^{\lbrack 0\rbrack} \right\rangle_{\theta}} \right)}} & {{Expression}\mspace{14mu} 8} \end{matrix}$

Here, x_(j) ^([0]) on the right side of Expression 8 denotes an input value supplied as training data by the supplying section 120, and <x_(j) ^([0])>_(θ) on the right side denotes a probability computed by using Expression 2 (1≤j≤I). The bias parameter b_(j) for each input node (1≤j≤I) may be adjusted and updated as denoted by the following expression by using Expression 8. Note that a coefficient c is a parameter predetermined by the user or the like.

$\begin{matrix} \left. b_{j}\leftarrow{b_{j} + {c\frac{1}{\tau}\left( {x_{j}^{\lbrack 0\rbrack} - \left\langle X_{j}^{\lbrack 0\rbrack} \right\rangle_{\theta}} \right)}} \right. & {{Expression}\mspace{14mu} 9} \end{matrix}$

That is, the learning processing section 150 adjusts the bias parameter b_(j) so as to increase the conditional probability of the input value x_(j) ^([0]) of the node of the 0-th input layer occurring, on a condition that the history x^((−T, −1]) of the common layer 12 has occurred. The learning processing section 150 may iteratively perform updating of the bias parameter b_(j) denoted by Expression 9 and computing of the probability <x_(j) ^([0])>_(θ) denoted by Expression 2, to determine the bias parameter b_(j). The learning processing section 150 stops updating the bias parameter b_(j) and determines the bias parameter b_(j) if a difference in the bias parameter b_(j) before and after updating is less than or equal to a predetermined threshold. If a bias parameter b_(j) is also set for a hidden node, the learning processing section 150 may determine the bias parameter b_(j) of the hidden node in the same manner.

Alternatively, the learning processing section 150 may decide upon the bias parameter b_(j) by iteratively updating the bias parameter b_(j) a predetermined number of times. If a difference in the bias parameter b_(j) before and after updating is greater than or equal to the predetermined threshold even after the bias parameter b_(j) has been updated the predetermined number of times, the learning processing section 150 may stop updating the bias parameter b_(j) and inform the user that the parameter does not converge.

Similarly, when updating the learning parameter u_(i,j,k), the learning processing section 150 may determine the direction of a change in the learning parameter u_(i,j,k) by using the following expression.

$\begin{matrix} {{\frac{\partial}{\partial u_{i,j,k}}\log\;{P_{\theta}\left( x_{\lbrack{1,I}\rbrack}^{\lbrack 0\rbrack} \middle| x^{({{- T},{- 1}})} \right)}} = {\tau^{- 1}{\alpha_{i,j,k}\left( {x_{j}^{\lbrack 0\rbrack} - \left\langle X_{j}^{\lbrack 0\rbrack} \right\rangle_{\theta}} \right)}}} & {{Expression}\mspace{14mu} 10} \end{matrix}$

In addition, when updating the learning parameter v_(i,j,l), the learning processing section 150 may determine the direction of a change in the learning parameter v_(i,j,l) by using the following expression.

$\begin{matrix} {{\frac{\partial}{\partial v_{i,j,l}}\log\;{P_{\theta}\left( x_{\lbrack{1,I}\rbrack}^{\lbrack 0\rbrack} \middle| x^{({{- T},{- 1}})} \right)}} = {{{- \tau^{- 1}}{\beta_{i,j,l}\left( {x_{j}^{\lbrack 0\rbrack} - \left\langle X_{j}^{\lbrack 0\rbrack} \right\rangle_{\theta}} \right)}} - {\tau^{- 1}{\gamma_{j,l}\left( {x_{i}^{\lbrack 0\rbrack} - \left\langle X_{i}^{\lbrack 0\rbrack} \right\rangle_{\theta}} \right)}}}} & {{Expression}\mspace{14mu} 11} \end{matrix}$

In the same manner as the updating of the bias parameter b_(j), the learning processing section 150 may iteratively perform updating of the learning parameters u_(i,j,k) and v_(i,j,l) corresponding to the input nodes (1≤j≤I) and computing of the probability <x_(j) ^([0])>_(θ) to determine the learning parameters u_(i,j,k) and v_(i,j,l) corresponding to the input nodes (1≤j≤I). Alternatively, the learning processing section 150 may iteratively perform an operation for updating the parameter set θ=(b_(j), u_(i,j,k), v_(i,j,l)) and then computing the probability <x_(j) ^([0])>_(θ) denoted by Expression 2 to determine the parameter set θ=(b_(j), u_(i,j,k), v_(i,j,l)).

As described above, the learning processing section 150 according to the present embodiment can decide upon the learning parameters u_(i,j,k) and v_(i,j,l) and the bias parameter b_(j) through learning. The learning apparatus 100 may then determine whether to continue learning (S360). The learning apparatus 100 may continue learning until it performs the learning process a predetermined number of times, or may continue learning until a stop command is input by the user. Alternatively, the learning apparatus 100 may continue learning until it can no longer acquire time-series data.

If the learning apparatus 100 continues learning (e.g., S360: YES), the process may return to step S310, in which the acquiring section 110 acquires the next time-series data, and the learning apparatus 100 may then perform learning of the model 10 based on the next time-series data. For example, the supplying section 120 supplies the 0-th input layer with the next image data in the image data acquired by the acquiring section 110. Furthermore, the storage section 130 samples the values of the hidden layers and stores these values in the 0-th hidden layer. Then, values held in the t-th common layer (−T<t<0) may be supplied to the (t−1)-th common layer. The values held in the (−T+1)-th layer may be deleted. The learning apparatus 100 may perform learning by using image data supplied to the layers from the 0-th input layer to the (−T+1)-th input layer as training data and using the values stored in the layers from the 0-th hidden layer to the (−T+1)-th hidden layer.

In this way, the supplying section 120 and the storage section 130 may sequentially acquire new input values x_(j) ^([0]) at the next time point corresponding to the nodes of the 0-th common layer in the model 10. Then, the computing section 140 may compute a conditional probability <x_(j) ^([0])>_(θ) of the new input value x_(j) ^([0]) on a condition that the history has occurred for each common layer before the next time point. The learning processing section 150 may adjust the weight parameter so as to increase the conditional probability of the new input value occurring on the condition that this history has occurred.

If the learning processing section 150 stops learning (e.g., S360: NO), the learning processing section 150 may output the learning parameters u_(i,j,k) and v_(i,j,l) and the bias parameter b_(j) that have been determined and store the parameters in the external database 1000 or the like.

As described above, the learning apparatus 100 according to an embodiment may be configured to apply, to time-series input data that is input in time series, a model having a total of T layers by associating one time point with the 0-th common layer and an input data sequence before the one time point with T−1 layers. The learning apparatus 100 may be configured to apply a model having hidden nodes to each common layer 12. That is, the learning apparatus 100 may be configured to form a time-evolution Boltzmann machine that predicts input data at one time point on the basis of the input data sequence and hidden node values.

The learning apparatus 100 is able to learn the model by computing a conditional probability of the input value x_(j) ^([0]) at the one time point occurring, based on the input value x^((−T, −1]), which is a history, for a model that takes time evolution into consideration. Furthermore, since the learning apparatus 100 learns the model using hidden nodes in addition to the time-series input data, the expressive ability and learning ability can be improved.

A description has been given of the learning apparatus 100 according to an embodiment that sequentially acquires new input data from time-series input data and adjusts the weight parameter for each input data acquired. Instead of this configuration, the learning apparatus 100 may acquire time-series input data of a predetermined duration and then adjust the weight parameters. For example, the learning processing section 150 adjusts the weight parameters collectively for a plurality of time points in response to acquisition of new input data at a plurality of time points corresponding to D layers.

FIG. 4 shows an example of structures of time-series data and training data for use in learning in an embodiment. In FIG. 4, the horizontal axis denotes time. FIG. 4 shows an example in which the learning apparatus 100 uses time-series data y^([1, L]) having a duration L that is longer than a duration T of time-series data y^([1, T]) used as training data by the learning apparatus 100 during learning. In this case, the learning processing section 150 may be configured to adjust weight parameters for a plurality of time points all together, in response to input data at a plurality of time points being newly acquired.

The learning apparatus 100 first performs learning using, as first training data, a time-series data segment of the time-series data from a time 1 to a time T. In this case, as described in FIG. 3, the learning apparatus 100 may perform learning by setting the time-series data and corresponding hidden layer values from the time 1 to the time T as each input value of the common layer 12 at the one time point in order, and incrementally shifting the time points one at a time toward the future. The learning apparatus 100 may use data at a time T as each input value x_(j) ^([0]) at the one time point, and continue learning until the time-series data from the time 1 to a time T−1 becomes the input data sequence x^((−T, 1]) (e.g., the history).

Next, the learning apparatus 100 performs learning using, as second training data, a time-series data segment of the time-series data from a time 2 to a time T+1. The learning apparatus 100 may sequentially use each of D pieces of data in the second training data as the input value x_(j) ^([0]) at the one time point. In this case, the learning apparatus 100 may shift the time point in the interval from the time 2 to the time T+1 one time point at a time toward the future and use, as the history, the corresponding time-series data and hidden nodes of the interval from the time 2 to the time T. In this way, the learning apparatus 100 may adjust the parameters D times for the D input values x_(j) ^([0]) and the corresponding D histories. That is, the learning apparatus 100 may use a stochastic gradient technique in which the learning method described with Expressions 8 to 11 is performed.

Alternatively, the learning apparatus 100 may acquire D time-series data sets, generate a plurality of training data sets from time-sequence data segments of a duration of L, and collectively perform learning for D layers. Specifically, the learning apparatus 100 may perform the stochastic gradient technique described using Expressions 8 to 11 collectively for D layers, by using the following expression.

$\begin{matrix} \left. \theta\leftarrow{\theta + {\eta{\sum\limits_{x \in {\{{y^{({t,{t + T - 1}}\rbrack}|{0 \leqq t < D}}\}}}^{\;}{{\nabla_{\theta}\log}\;{P_{\theta}\left( x_{\lbrack{1,I}\rbrack}^{\lbrack 0\rbrack} \middle| x^{({{- T},{- 1}})} \right)}}}}} \right. & {{Expression}\mspace{14mu} 12} \end{matrix}$

FIG. 5 shows a first modification of the learning apparatus 100 according to an embodiment. Components of the learning apparatus 100 shown in FIG. 5 that perform substantially the same operations as those of the learning apparatus 100 according to the embodiment illustrated in FIG. 1 are denoted by the same reference numerals, and a description thereof is omitted. In a case where time-series data of a duration L, such as described in FIG. 4, is provided, the learning apparatus 100 according to the present modification may be configured to efficiently update parameters by using FIFO memories and learn a model corresponding to the time-series input data. The learning apparatus 100 according to the present modification further includes FIFO memories 160 and an updating section 170.

Each of the FIFO memories 160 may sequentially store input data and output the stored data after a predetermined number of storages have been performed. Each of the FIFO memories 160 may be a memory that first outputs data that has been stored first (e.g., FIFO: First In, First Out).

Each of the FIFO memories 160 may sequentially store an input value of the common layer 12 and output the input value after a predetermined number of storages have been performed. The learning apparatus 100 may include a plurality of FIFO memories 160, the number of which is greater than or equal to the number of nodes n of the model. The plurality of FIFO memories 160 is desirably provided to have a one-to-one correspondence with the plurality of nodes of the common layer 12. That is, each of the plurality of FIFO memories 160 may be provided in a manner to store a history for a respective node of the common layer 12 or to update the history thereof.

The plurality of FIFO memories 160 are connected to the acquiring section 110 and the storage section 130, and sequentially store input values corresponding to new input data of the common layer 12. The plurality of FIFO memories 160 are also connected to the updating section 170 and sequentially supply the data stored therein to the updating section 170.

The updating section 170 may be configured to update a plurality of update parameters that are based on the hidden nodes and the input data sequence of the time-series input data before the one time point, from values at a prior time point to values at the one time point, on the basis of values of the update parameters and values of the hidden nodes and input values corresponding to the input data to be reflected next. The updating section 170 may update the update parameters by using values input to the FIFO memories 160 and values output from the FIFO memories 160. The updating section 170 may be connected to the acquiring section 110 and the storage section 130, and may receive values input to the FIFO memories 160. Alternatively, the updating section 170 may receive values input to the FIFO memories 160 from the acquiring section 110 via the supplying section 120.

Here, the update parameters are α_(i,j,k) and γ_(i,l) shown in Expressions 5 and 7. In this case, the update parameters are based on input values i (1≤i≤I) corresponding to input data of the input data sequence at each time point and the predefined parameters λ_(k) ^(t1) and μ_(l) ^(t2) of the weight parameter W_(ij) ^([δ]) between this input value i and the target input node j (1≤j≤I) or hidden node j (I+1≤j≤I+H), for example. As another example, the update parameters are based on the hidden node i (I+1≤i≤I+H) at each time point and the predefined parameters λ_(k) ^(t1) and μ_(l) ^(t2) of the weight parameter W_(ij) ^([δ]) between this hidden node i and the target input node j (1≤j≤I) or hidden node j (I+1≤j≤I+H), for example.

The update parameters may be updated every time the acquisition of the time-series input data by the acquiring section 110 and the storage of the sampling values by the storage section 130 are performed sequentially. The above-described learning apparatus 100 according to the present modification may be configured to learn a modification of the model 10. The modification of the model 10 is described with reference to FIG. 6.

FIG. 6 shows a modification of the model 10 according to an embodiment. The model 10 according to the modification needs not have the layered structure including T layers shown in FIG. 2. FIG. 6 shows an example of a model 10 corresponding to one of the FIFO memories 160. Accordingly, the overall configuration of the model 10 according to the present modification includes a storage area that is equivalent to the 0-th common layer in FIG. 2 including the training data, and a number of the configurations illustrated in FIG. 6 equal to the number of nodes n (=I+H). Neurons i and j and a FIFO sequence 20 of the model 10 according to the present modification are described below.

The neuron i may be equivalent to the input terminal of the FIFO memory 160. An input value y_(i) ^([t]) (1≤i≤I) of each node in the input data of the input data sequence at each time point t and a corresponding value y_(i) ^([t]) among the values y_(i) ^([t]) (I+1≤i≤I+H) of the hidden nodes at each time point are sequentially input to the neuron i. The neuron i may set the value y_(i) ^([t]) input thereto as the current input value. Then, at a time point t+1, the neuron i may supply the input value y_(i) ^([t]) input at the time point t to the updating section 170 and to the FIFO sequence 20 as the previous input value and may hold the input value y_(i) ^([t+1]) at the time point t+1 as the current input value.

The FIFO sequence 20 may store d_(ij)−1 of the latest input values received from the neuron i. The FIFO sequence 20 may supply the d_(ij)−1 input values stored therein to the updating section 170. The updating section 170 may be configured to compute the values of the update parameters denoted by Expression 6 by using the input values supplied by the FIFO sequence 20. If the FIFO sequence 20 holds input values from the time point t−1 to the time point t−d_(ij)+1, the FIFO sequence 20 is denoted by the following expression. q _(i,j)≡(y _(i) ^([t−1]) ,y _(i) ^([t−d) ^(ij) ^(+2]) ,y _(i) ^([t−d) ^(ij) ^(+1]))  Expression 13:

After the input value y_(i) ^([t1]) is input to the neuron i at the time point t1, the FIFO sequence 20 may store the input value y^(i[t1]) up until a time point t3 (=t1+d_(ij)−1) which is a predetermined time period d_(ij)−1 after the next time point t2 (=t1+1) of the time point t1. At the next time point t4 (=t3+1=t1+d_(ij)), the FIFO sequence 20 may supply the input value y_(i) ^([t1]) to the neuron j. The input value y_(i) ^([t1]) supplied to the neuron j at the time point t4 is immediately supplied to the updating section 170 at the time point t4. However, the input value y_(i) ^([t1]) that the neuron j has received from the FIFO sequence 20 at the time point t4 does not serve as an input for the neuron j, and the input value y_(i) ^([t4]) may be input to the neuron j at the time point t4.

The neuron j may be equivalent to the output terminal of the FIFO memory 160, and the neuron j may receive the input value y_(i) ^([t1]) input to the neuron i at the time point t1, via the FIFO sequence 20 after the time period d_(ij), e.g., at the time point t1+d_(ij). That is, the model 10 from the neuron i to the neuron j via the FIFO sequence 20 may correspond to the FIFO memory 160 that stores d_(ij) pieces of input data. In addition, the neuron i of the model 10 according to the modification may correspond to, for example, a node for an input data sequence such as a node i of the (−δ)-th common layer of the model 10 shown in FIG. 2, and in this case the neuron j may correspond to, for example, the node j of the 0-th common layer. At the time point t1+d_(i,j), the neuron j may supply the received input value y_(i) ^([t1]) to the updating section 170.

As described above, the model 10 according to the present modification may supply the input values at the time point t−1 and the time point t-d_(ij)+1 to the updating section 170 at the time point t. In this way, the updating section 170 can update the update parameters by adding the corresponding input value in the input data to be reflected next to the update parameters for the time point before the one time point, and then multiplying the resulting sum by a predetermined constant. Note that the update parameters denoted by Expression 8 may be computed in accordance with Expression 8 by using the input values stored in the FIFO sequence 20 that are supplied to the updating section 170.

For example, the update parameter γ_(i,l) denoted by Expression 7 can be updated by using the input values supplied to the updating section 170 and the second predefined parameter. Specifically, the updating section 170 can compute the update parameter γ_(i,l) to be used in the current learning by performing computing at the time point t according to the following expression by using the prior update parameter γ_(i,l) and the input value γ_(i) ^([t−1]) received from the neuron i at the time point t. γ_(i,l)←μ_(l)(γ_(i,l) +y _(i) ^([t−1]))  Expression 14:

FIG. 7 shows an example of a temporal change in the update parameter γ_(i,l) according to the present embodiment. FIG. 7 shows an example in which values greater than 0 (for example, 1) are input to the neuron i as the input value at time points t−5, t−2, and t−1, and these input values are supplied to the updating section 170 at time points t−4, t−1, and t. The second predefined parameter μ_(l) is a parameter whose value gradually decreases as the time point difference increases. Accordingly, the update parameter γ_(i,l) computed by the updating section 170 tends to decrease as time passes from when the input value of 1 is input to when the next input is given.

The update parameter α_(i,j,k) denoted by Expression 5 can be updated by using the input values supplied to the updating section 170 and the first predefined parameter λ_(k). Specifically, the updating section 170 can compute the update parameter α_(i,j,k) to be used in the current learning by performing computing at the time point t according to the following expression, by using the prior update parameter α_(i,j,k) and the input value γ_(i) ^([t−dij]) received from the neuron j at the time point t. α_(i,j,k)←λ_(k)(α_(i,j,k) +y _(i) ^([t−d) ^(ij) ^(]))  Expression 15:

FIG. 8 shows an example of a temporal change in the update parameter α_(i,j,k) according to the present embodiment. FIG. 8 shows an example in which values greater than 0 (for example, 1) are supplied to the neuron j as the input value at time points t−3, t−1, and t. The first predefined parameter λ_(k) is a parameter whose value gradually decreases as the time point difference increases. Accordingly, the update parameter α_(i,j,k) computed by the updating section 170 tends to decrease as time passes from when the input value of 1 is input to when the next input is given.

As described above, the learning apparatus 100 according to the present modification can update the update parameters α_(i,j,k) and γ_(i,l) by applying the model 10 shown in FIG. 6 using the FIFO memories 160 and the updating section 170. Note that the updating section 170 can apply the model 10 according to the present modification, for example, by acquiring the input values x_(i) ^([t−1]) at the time point t−1 from the input data input to the FIFO memories 160 and acquiring the input values x_(i) ^([t−dij]) at the time point t−d_(ij) from the output of the FIFO memories 160.

In addition, the learning apparatus 100 may update the parameter β_(i,j,l) through substantially the same operation as the operation described in FIG. 3. Specifically, the computing section 140 can compute the parameter β_(i,j,l) by determining the sum of products of the second predefined parameter μ_(l) and the input value x_(i) (y_(i) in the present modification) for time points from t−1 to t−d_(ij)+1 as indicated by Expression 6.

In this way, the computing section 140, according to the present modification, can compute, by using the plurality of update parameters, conditional probabilities of input data values at one time point on the condition that the hidden node values and input data sequence have occurred. Then, the learning processing section 150 can determine the learning parameters u_(i,j,k) and v_(i,j,l) and the bias parameter b_(j) by performing substantially the same operation as the operation described in FIG. 3.

In other words, the learning apparatus 100, according to the present embodiment, can determine the weight parameter and bias parameters in a manner to increase the probability of predicting the input value to be input to the input layer 14 of the common layer 12, based on the past values that have been input to the common layer 12 of the model 10 before the one time point. Furthermore, the learning apparatus 100 can improve the prediction accuracy, the expressive ability, the learning efficiency, and the like of the input values input to the input layer 14 by having the common layer 12 include the hidden layer 16 in addition to the input layer 14.

The learning apparatus 100, according to the present embodiment described above, is an example in which a value that is unrelated to the prediction made by the learning apparatus 100 is sampled and input as the hidden node value to be input to the hidden layer 16. Instead, the learning apparatus 100 may determine the hidden node value by using a history of the conditional probability of the values of the nodes of the common layer 12. The learning apparatus 100 may determine the weight parameter to a hidden node using this conditional probability history. The learning apparatus 100 can improve the prediction accuracy by using the conditional probability history of nodes of the common layer 12 to determine the weight parameter to the hidden node and the hidden node value.

In this case, the computing section 140 may compute the conditional probability p_(j,t) of the value of a node j of the common layer 12 at one time point t based on the values input to the corresponding node j of the common layer 12 at each time point before the one time point t, and store this conditional probability in the storage section or the like. In addition to the computation of the conditional probability of each input value of the input layer 14 at the one time point described above, the computing section 140 may compute the conditional probability of each hidden node in the layer 16 at the one time point in the same manner. That is, the computing section 140 may use the plurality of update parameters to compute the conditional probability of the value of each hidden node and each input data value at the one time point on the condition that an input data sequence has occurred. Here, the computing section 140 may store the conditional probability p_(j,t) in a FIFO or the like.

The computing section 140 may be configured to compute a total likelihood, after the learning by the learning apparatus 100 has continued. The computing section 140 computes the total likelihood p_(j) as shown in the following expression, based on the conditional probabilities p_(j,t−K+1), p_(j,t−K+2), . . . , p_(j,t) computed by K instances of learning from the time point t−K+1 to the time point t, for example. The total likelihood p_(j) in Expression 16 indicates a total sum of the conditional probabilities, as an example, but the total likelihood p_(j) may be at least one of a sum, weighted sum, product, or weighted product of the conditional probabilities. Furthermore, K may be an integer greater than or equal to 2, and if the computing section 140 stores the conditional probabilities p_(j,t) in a FIFO or the like, the length of the FIFO sequence may be equal to the value of K.

$\begin{matrix} {p = {\sum\limits_{s = {t - k + 1}}^{t}p_{s}}} & {{Expression}\mspace{14mu} 16} \end{matrix}$

The computing section 140 may supply the total likelihood p_(j) to the storage section 130. The storage section 130 may sample the values x_(j) ^([t]) of the hidden nodes of the hidden layer 16 at the one time point, based on the most recent likelihood p_(j,t). That is, the storage section 130 according to the present embodiment may be configured to sample the value of each hidden node at the one time point, by using the conditional probability of the value of each hidden node at the one time point. For example, the storage section 130 samples the values of the hidden nodes based on the history of the conditional probabilities computed by the computing section 140. That is, the storage section 130 may sample the values of the hidden nodes after the learning operation of the learning processing section 150 has been performed a plurality of times. The storage section 130 may store a value of 0 in the hidden nodes as the sampling value, until the learning operation of the learning processing section 150 has been performed a plurality of times.

The storage section 130 may store a value of 1 or 0 in the hidden node j as the sampling result, according to the result of a comparison between the value of the total likelihood p_(j) and a threshold value. In this way, when predicting the time series data to be input to the input layer 14, the storage section 130 can store a more preferable value as the hidden node value by performing sampling based on the history of past conditional probabilities.

The learning processing section 150 may be configured to determine the weight parameter based on the total likelihood p_(j). In this case, the learning processing section 150 may compute update amounts Δu_(i,j,k) and Δv_(i,j,k) for the learning parameters u_(i,j,k) and v_(i,j,k) in the weight parameter for one hidden node j at the one time point. For example, the learning processing section 150 may compute these update amounts Δu_(i,j,k) ^([t]) and Δv_(i,j,k) ^([t]) as shown in the following expression, based on the value x_(j) ^([t]) of the one hidden node j at the one time point t and on the conditional probability <X_(j) ^([t])> of the value of this hidden node j at the one time point t on the condition that the input data sequence has occurred (I+1≤j≤I+H). Δu _(i,j,k) ^([t])=α_(i,j,k) ^([t−1])(x _(j) ^([t]) −

X _(j) ^([t])

)  Expression 17: Δv _(i,j,l) ^((1)[t])=β_(i,j,l) ^([t−1])(

X _(j) ^([t])

−

x _(j) ^([t])) Δv _(i,j,l) ^((2)[t])=γ_(j,l) ^([t−1])(

X _(i) ^([t])

−

x _(i) ^([t]))

Here, the update amount Δv_(i,j,k) ^([t]) is equal to Δu_(i,j,k) ^((1)[t])+Δu_(i,j,k) ^((2)[t]). The conditional probability <X_(j) ^([t])> of the value of the hidden node j may be computed by the computing section 140 using Expression 2. The learning processing section 150 may store the computed update amounts Δv_(i,j,k) ^([t]), Δu_(i,j,k) ^((1)[t]), and Δu_(i,j,k) ^((2)[t]) in the storage section or the like. The learning processing section 150 may be configured to store the update amounts Δv_(i,j,k) ^([t]), Δu_(i,j,k) ^((1)[t]), and Δu_(i,j,k) ^((2)[t]) computed for one time point in the FIFO sequence. That is, the learning processing section 150 may be configured to update the learning parameters based on update amounts computed in the past.

The learning processing section 150 changes the ratio by which the update amounts are reflected in the learning parameters, according to the conditional probability of input data occurring at a following time point that is after the one time point t, for example. In this case, the learning processing section 150 may change the ratio by which the update amounts are reflected in the learning parameters according to the conditional probability of a hidden node value occurring at a plurality of following time points that are after the one time point. In the present embodiment, an example is described in which the learning processing section 150 changes the ratio by which the update amounts are reflected in the learning parameters according to the total likelihood p_(j) after the one time point.

The learning processing section 150 may update the learning parameters as shown in the following expression, based on the total likelihood p_(j) computed by K instances of learning from the time point t−K+1 to the time point t and on the update amount at the time point t−K+1, for example. Here, K may be an integer greater than or equal to 2, and if the learning processing section 150 stores the update amounts in the FIFO sequence, the length of the FIFO sequence may be equal to the value of K. u _(i,j,k) ←u _(i,j,k)+η₁ p _(j) Δu _(i,j,k) ^([t−k+1])  Expression 18: v _(i,j,l) ←v _(i,j,l)+η₁ p _(j)(Δv _(i,j,l) ^((1)[t−k+1]) +Δv _(i,j,l) ^((2)[t−k+1]))

Here, η₁ may be a constant for adjusting the update amount. Alternatively, η₁ may be a coefficient whose value becomes smaller according to an increase in the number of updates. Yet further, η₁ may have a value of substantially 1 at the stage when the learning processing section 150 begins learning, and may be a coefficient whose value becomes smaller according to the amount of learning occurring as time progresses from the time point t. For example, η₁=η₁₀/t². Furthermore, η₁ may be a coefficient whose value becomes smaller according to the update amount. For example, η₁=η₁₀/(ΣΔu_(i,j,k) ²)^(1/2). Here, η₁₀ may be a predetermined constant.

In the manner described above, the learning processing section 150 may update the learning parameters of a hidden node at one time point according to the conditional probabilities computed at time points before the one time point. In this way, the learning apparatus 100 can more strongly reflect the update amounts at time points before the one time in the learning parameters, in response to the predicted probability of an input value of an input node being large due to the weight parameters at time points before the one time point. That is, the learning apparatus 100 can update the weight parameters of the hidden nodes in a manner to increase the conditional probabilities.

If a FIFO sequence is used to perform an update of such a weight parameter, the learning processing section 150 may extract from the FIFO sequence the update amounts Δv_(i,j,k) ^([t−K+1]), Δu_(i,j,k) ^((1)[t−K+1]), and Δu_(i,j,k) ^((2)[t−K+1]) of a past time point (e.g., t−K+1) at the following time point (e.g., t) or a time point thereafter. The learning processing section 150 may update the ratio by which the update extracted from the FIFO sequence are reflected in the learning parameters according to the conditional probabilities of hidden node values occurring at the following time point t. For example, the learning processing section 150 may multiply the total likelihood p_(j) respectively by each update amount. In this way, the learning processing section 150 can efficiently perform the update of the weight parameters as described above.

The learning apparatus 100, according to the embodiment described above, is an example in which the input value x_(j) ^([0]) of each node is a binary value of 1 or 0. The learning apparatus 100 calculates the conditional probability for such binary data, but if time-series data made up of multi-values or real values is input as-is, the calculation result of the calculating section 140 cannot be handled as a probability. However, by using a value corresponding to a node value made up of a multi-value or a real value, the learning apparatus 100 can be configured to process multi-value or real number time-series data. The following describes such a learning apparatus 100, as a second modification of the learning apparatus 100.

The learning apparatus 100 of the second modification can learn a model corresponding to multi-value or real number time-series input data, using substantially the same configuration as shown in FIG. 1 or FIG. 5. In the present embodiment, the second modification of the learning apparatus 100 is described using the configuration shown in FIG. 5. In this case, the acquiring section 110 may be configured to acquire time-series input data that is a time series of input data including a plurality of input values x_(j) ^([0]). The input values x_(j) ^([0]) may include multi-values or real values. The operation of the acquiring section 110 may be the same as the operation of the acquiring section 110 already described above, except that the time-series input data that is acquired includes multi-values or real values, and therefore this operation is omitted from the description.

The supplying section 120 may be configured to supply a plurality of the nodes of the model with a plurality of input values corresponding to pieces of input data at one time point in the time-series input data. In other words, the supplying section 120 may be configured to supply a plurality of input values to the input side of the FIFO memory 160. The operation of the supplying section 120 may be the same as the operation of the supplying section 120 already described above, except that the input values being stored include multi-values or real values, and therefore this operation is omitted from the description. Similarly, the operation of the storage section 130 may be the same as the operation of the storage section 130 already described above, except that the input values being stored include multi-values or real values, and therefore this operation is omitted from the description. Furthermore, the operation of the FIFO memory 160 may be the same as the operation of the FIFO memory 160 already described above, except that the input values being stored include multi-values or real values, and therefore this operation is omitted from the description.

The updating section 170 may be configured to update the update parameter using the values input to the FIFO memory 160 and the values output from the FIFO memory 160. The operation of the updating section 170 may be the same as the operation of the updating section 170 already described above, except that the input data includes multi-values or real values, and therefore this operation is omitted from the description.

The calculating section 140 may be configured to calculate the node value of each of a plurality of nodes of the model, instead of calculating a conditional probability or in addition to calculating a conditional probability. The calculating section 140 may be configured to, if the input value of a node is binary, calculate the conditional probability corresponding to this node. The calculation of the conditional probability performed by the calculating section 140 may be the same as the operation of the calculating section 140 already described above, and therefore this operation is omitted from the description.

The calculating section 140 may be configured to, if the input value of a node is a multi-value or a real value, calculate the node value corresponding to this node. The calculating section 140 may be configured to calculate an average value of possible values of a node on the condition that the input data sequence has occurred, as the node value of this node. The calculating section 140 may calculate each node value corresponding to each node at one time point, based on the input data series before the one time point in the time-series input data and the weight parameter in the model.

The learning processing section 150 may be configured to adjust the weight parameters in the model. The learning processing section 150 may be configured to, if the input value of a node is a binary value, further increase the conditional probability of the input data of this node occurring at one time point on the condition that the input data sequence has occurred. The adjustment of the weight parameters performed by the learning processing section 150 may be the same as the operation of the learning processing section 150 already described above, and therefore this operation is omitted from the description.

The learning processing section 150 may be configured to, if the input value of a node is a multi-value or a real value, update the weight parameter using the corresponding input value and the calculated error of the node value at one time point. Furthermore, the learning processing section 150 may be configured to further update the variance parameter for indicating variance in the probability distribution of the input value using the corresponding input value and the calculated error of the node value at one time point.

As described above, the learning apparatus 100 of the second modification may be configured to, if the input value of a node is a multi-value or a real value, express the input value of this node using the variance and the average value of the possible values of this node. The following describes the operation of the learning apparatus 100 of the second modification in a case where the input value of a node is a multi-value or a real value.

FIG. 9 shows an operational flow of the learning apparatus 100 of a second modification according to the present embodiment. FIG. 9 shows an example in which the learning apparatus 100 of the second modification operates according to a multi-value or real value input value, using the model shown in FIG. 6. Specifically, the learning apparatus 100 may be configured to calculate each propagation value (α_(i,j,k), β_(i,j,l), and γ_(i,l)) weighted according to the passage of time points for each node value at a plurality of time points before the one time point. The calculating section 140 then propagates each propagation value to the plurality of nodes in the model in correspondence with the plurality of input values x_(j) ^([t]) at the one time point t. In the present embodiment, a node that is a propagation destination corresponding to real number time-series input data is referred to as a first node. Here, the propagation destination node may be a node corresponding to a neuron j in the model 10.

First, the acquiring section 110 may acquire the multi-value or real value time-series data (S410). The acquiring section 110 may acquire the time-series data of an interval corresponding to a time from one time point to a time point that is a predetermined first number of time points before the one time point in the model 10. For example, the acquiring section 110 acquires T pieces of real value data arranged in time series.

Next, the supplying section 120 may supply the input values corresponding to the input data at the one time point in the time-series input data respectively to the input nodes corresponding to the 0-th layer of the model 10 (S420). Here, the input data supplied to each input node at the one time point t is x_(j) ^([t]) (1<j<I). The supplying section 120 may supply the input values corresponding to the input data series at time points before the one time point to the FIFO sequence 20 of the model 10 in order from the oldest time point. For example, a real value history up to when the input data reaches the input data x_(j) ^([t]) in the time-series data is input to the FIFO sequence 20.

If hidden nodes are present in the model 10, the storage section 130 may sample the values of the hidden nodes corresponding to the one time point and respectively store the sampled values in the corresponding one or more hidden nodes j (I+1<j<I+H) (S430). The storage section 130 may sample multi-values or real number values in correspondence with the time-series input data.

The calculating section 140 may calculate the node value B_(j) ^([t]) corresponding to each input value x_(j) ^([t]) (1<j<I) of the input nodes at the one time point, based on the input values x_(j) ^([−T, −1]) of the plurality of nodes including the hidden nodes and the weight parameters (S440). The calculating section 140 may be configured to calculate the node values B_(j) ^([t]) of the first nodes among the plurality of nodes by using each propagation value propagated to a first node. The calculating section 140 may calculate the node value B_(j) ^([t]) as shown in the following expression, based on each propagation value and the weight parameters.

$\begin{matrix} {B_{j}^{\lbrack t\rbrack} = {b_{j} + {\sum\limits_{i = 1}^{I + H}\left( {{\sum\limits_{k \in K}^{\;}{u_{i,j,k}\alpha_{i,j,k}}} - {\sum\limits_{l \in L}^{\;}{v_{i,j,l}\beta_{i,j,l}}} - {\sum\limits_{l \in L}^{\;}{v_{j,i,l}\gamma_{i,l}}}} \right)}}} & {{Expression}\mspace{14mu} 19} \end{matrix}$

In this way, when the time-series input data includes real values, the calculating section 140 may calculate the node value B_(j) ^([t]) of a corresponding first node by using Expression 19, which is a portion of Expression 3. The node value B_(j) ^([t]) of a first node at the one time point t calculated by the calculating section 140 is an average value of the possible values of this first node at the one time point t.

Next, the learning processing section 150 may update the weight parameter used when calculating each propagation value propagated to a first node by using the error (e.g., x_(j) ^([t])−B_(j) ^([t])) between the corresponding input value x_(j) ^([t]) and the calculated node value B_(j) ^([t]) at the one time (S450). Here, if the learning is performed using the gradient technique corresponding to the stochastic gradient technique described in Expressions 8 to 11, the learning processing section 150 may update the weight parameter based on a normal distribution. In other words, the learning processing section 150 may be configured to update the weight parameter by using the variance parameter σ for indicating the variance in the probability distribution p(x_(j) ^([t])) of the input values x_(j) ^([t]) of the first nodes. In this case, the probability distribution p(x_(j) ^([t])) of the input values x_(j) ^([t]) is expressed as shown below.

$\begin{matrix} {{p\left( x_{j}^{\lbrack t\rbrack} \right)} = {\frac{1}{\sqrt{2\pi}\sigma}\exp\left\{ {- \frac{\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)^{2}}{2\sigma^{2}}} \right\}}} & {{Expression}\mspace{14mu} 20} \end{matrix}$

The learning processing section 150 may update the learning parameters as shown in the following expression. As shown in the following expression, when updating the weight parameters, the learning processing section 150 may set the update amount of the weight parameters to be smaller when the variance parameter σ is larger. Furthermore, in the same manner as η1, η2 may be a constant for adjusting the update amount or may be a coefficient whose value becomes smaller according to an increase in the number of updates.

$\begin{matrix} {\mspace{79mu}{\left. u_{i,j,k}\leftarrow{u_{i,j,k} + {\eta_{2}\frac{1}{\sigma^{2}}{\alpha_{i,j,k}^{\lbrack{t - 1}\rbrack}\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)}}} \right.\left. v_{i,j,l}\leftarrow{v_{i,j,l} + {\eta_{2}\frac{1}{\sigma^{2}}{\beta_{i,j,l}^{\lbrack{t - 1}\rbrack}\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)}} + {\eta_{2}\frac{1}{\sigma^{2}}{\gamma_{i,l}^{\lbrack{t - 1}\rbrack}\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)}}} \right.}} & {{Expression}\mspace{14mu} 21} \end{matrix}$

The learning processing section 150 may be configured to further update the variance parameter σ. For example, the learning processing section 150 may update the variance parameter σ by using the error between the corresponding input value and the calculated node value at the one time for a first node. The learning processing section 150 may be configured to, when updating the variance parameter σ, update the variance parameter σ based on the mathematical square of the error for each of a plurality of nodes, as shown in the expression below. In other words, the variance parameter σ may be a parameter that is common to a plurality of input values.

$\begin{matrix} \left. \sigma\leftarrow{\sigma + {\eta_{2}\left( {{- \frac{I}{\sigma}} + {\frac{1}{\sigma^{3}}{\sum\limits_{j = 1}^{I}\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)^{2}}}} \right)}} \right. & {{Expression}\mspace{14mu} 22} \end{matrix}$

Next, the updating section 170 may update the update parameter (S460). For example, the updating section 170 may update the update parameter as shown in the expression below. The updating section 170 may be configured to update the update parameter α_(i,j,k) based on the input value x_(j) ^([t5]) that is a first number d_(ij) of time points before the one time point (t5=t−d_(ij)). Here, the update parameter α_(i,j,k) is set as a first update parameter. The updating section 170 may be configured to, when updating the first update parameter α_(i,j,k), update the first update parameter α_(i,j,k) by multiplying a first coefficient by the sum of the original update parameter α_(i,j,k) ^([t−1]) and the input value x_(j) ^([t5]) at a time point that is the first number d_(ij) of time points before the one time point t. Here, the first predefined parameter λ_(k) is set as the first coefficient.

$\begin{matrix} \left. \left. \alpha_{i,j,k}^{\lbrack t\rbrack}\leftarrow{\lambda_{k}\left( {\alpha_{i,j,k}^{\lbrack{t - 1}\rbrack} - x_{i}^{\lbrack{t - d_{ij}}\rbrack}} \right)} \right.\left. \beta_{i,j,l}^{\lbrack t\rbrack}\leftarrow{\sum\limits_{s = {{- d_{ij}} + 1}}^{- 1}{x_{i}^{\lbrack{t + s}\rbrack}\mu_{l}^{s}}} \right.\gamma_{i,l}^{\lbrack t\rbrack}\leftarrow{\mu_{l}\left( {\gamma_{i,l}^{\lbrack{t - 1}\rbrack} + x_{i}^{\lbrack t\rbrack}} \right)} \right. & {{Expression}\mspace{14mu} 23} \end{matrix}$

Here, since the time-series input data is a real number, the update parameter is also a real value. For example, it is possible for the change over time of the update parameters α_(i,j,k) and γ_(i,l) to result in the input value being a negative value. In such a case, if the value of the input parameter at the time point when the input value became this negative value is less than an absolute value of this negative input value, the update parameter becomes a negative value. The update parameter calculated by the updating section 170 then exhibits a trend of increasing over time to draw near a value of 0.

In the manner described above, the learning processing section 150 can determine the learning parameters u_(i,j,k) and v_(i,j,l) and the variance parameter σ by learning from the real number input data time series. Furthermore, the updating section 170 can update the update parameters. The learning apparatus 100 may determine whether to continue this learning (S470). The learning apparatus 100 may continue learning until reaching a predetermined number of learning processes, or may instead continue learning until a stop command is input from the user. As another example, the learning apparatus 100 may continue learning until there is no more time-series data that can be acquired.

If learning continues (e.g., S470: Yes), the learning apparatus 100 returns the processing to step S410 and, if there is no more time-series data to be supplied to the FIFO sequence 20, the acquiring section 110 may acquire the next piece of time-series data and the learning apparatus 100 may learn the model 10 based on this next piece of time-series data. The supplying section 120 supplies the next piece of real number data in the time-series data acquired by the acquiring section 110 to the corresponding FIFO sequence 20, for example. The supplying section 120 may supply the FIFO sequence 20 with the data from the oldest time point in in the time-series data to be supplied to the FIFO sequence 20. The storage section 130 may sample the values of a hidden layer and supply these values to the corresponding FIFO sequence 20.

The calculating section 140 may calculate the node value B_(j) ^([t]) of the first node based on a value obtained by weighting the updated update parameter with the updated weight parameter. For example, the calculating section 140 may calculate the node value B_(j) ^([t]) by using the value u_(i,j,k)·α_(i,j,k) obtained by weighting the update first update parameter α_(i,j,k) with the first weight u_(i,j,k) included in the updated weight parameter. The learning processing section 150 and the updating section 170 may update each parameter based on the new node value B_(j) ^([t]).

In this way, the supplying section 120 and the storage section 130 may sequentially input new real values at subsequent time points to the FIFO sequences of the model 10. The calculating section 140 may calculate the new node value B_(j) ^([t]) on the condition that there is a history of real values before the subsequent time point. The learning processing section 150 and the updating section 170 may sequentially perform learning of the real number time-series input data by updating each parameter.

If the learning is ended (e.g., S470: No), the learning processing section 150 may output each determined parameter and store these parameters in the external database 1000 or the like. As described above, the learning apparatus 100 of the second modification may be configured to form a Boltzmann machine that considers time expansion for predicting input data at one time point by using the average value and variance for real number time-series input data.

The above describes an example in which the learning apparatus 100 of the second modification operates according to real number or multi-value time-series input data. In addition to this, the learning apparatus 100 may be configured to operate according to time-series input data including real numbers, multi-values, and binary values. In other words, at least a portion of the plurality of input values may be multi-values or real values.

For example, the following describes a case in which one portion of the input values consists of multi-values or real values and another portion of the input values consists of binary values. In this case, the learning apparatus 100 may calculate, for a first node associated with an input value that is a multi-value or real value among the plurality of input values, the node value B_(j) ^([t]) of this first node based on the operational flow shown in FIG. 9.

Furthermore, the learning apparatus 100 may calculate, for a second node associated with an input value that is a binary value among the plurality of input values, the conditional probability described in FIGS. 1 to 8. In other words, the learning apparatus 100 may be configured to calculate the conditional probability of the second node being one of the two values of the binary value based on each propagation value (α_(i,j,k), β_(i,j,l), and γ_(i,l)) propagated to the second node and the input value x_(j) ^([t]) at the one time point corresponding to the second node.

In other words, the learning apparatus 100 may be configured to calculate the conditional probability of the value of a second node at one time point on a condition that an input data time series has occurred, based on the input data time series before the one time point in the binary time-series input data and this weight parameter in the model. The learning processing section 150 may adjust the weight parameter in the model to increase the conditional probability of input data occurring at the one time point on a condition that the input data time series has occurred.

The calculating section 140 may be configured to calculate, for a second node associated with an input value that is a binary value among the plurality of input values, the expected value <X_(j) ^([t])> of the node value of the second node, based on each propagation value (α_(i,j,k), β_(i,j,l), and γ_(i,l)) propagated to the second node and the input value x_(j) ^([t]) at the one time point corresponding to the second node. In this case, the calculating section 140 may be configured to update the weight parameter used to calculate each propagation value propagated to the second node by using the error (x_(j) ^([t])−<X_(j) ^([t])>) between the corresponding input value x_(j) ^([t]) and the expected value <X_(j) ^([t])> of the node value of the second node at the one time point t. The update of the weight parameter is the same as the operation described in Expression 17, Expression 18, and the like, and therefore is omitted from the description.

In the manner described above, the learning apparatus 100 may perform a different process for each model by respectively using a FIFO sequence 20 corresponding to real number or multi-value time-series input data and a FIFO sequence 20 corresponding to binary time-series input data. In this way, the learning apparatus 100 can be used for time-series input data including real numbers, multi-values, and binary values, and can therefore increase the expressive power and learning power and can be applied in various fields, such as moving images, languages, and music.

The learning apparatus 100 of the second modification described above is an example in which the calculating section 140 calculates the node value B_(j) ^([t]) using Expression 19. However, the expression used by the calculating section 140 to calculate the node value B_(j) ^([t]) is not limited to this. For example, the calculating section 140 may calculate the node value B_(j) ^([t]) as shown by the following expression.

$\begin{matrix} {B_{j}^{\lbrack t\rbrack} = {b_{j} + {\sum\limits_{i = 1}^{I + H}\left( {{\sum\limits_{k \in K}^{\;}{u_{i,j,k}\alpha_{i,j,k}}} + {\sum\limits_{\delta = 1}^{d_{ij} - 1}{w_{ij}^{\lbrack\delta\rbrack}x_{i}^{\lbrack{t - \delta}\rbrack}}} - {\sum\limits_{l \in L}^{\;}{v_{j,i,l}\gamma_{i,l}}}} \right)}}} & {{Expression}\mspace{14mu} 24} \end{matrix}$

In this way, the calculating section 140 may calculate the node value B_(j) ^([t]) using a weight parameter including second weights w_(ij) corresponding to input values at two or more time points between the one time point t and a time point that is a first number d_(ij) of time points before the one time point t. Specifically, the calculating section 140 may calculate the node value B_(j) ^([t]) using the d_(ij)−1 input values x_(i) ^([t−δ]) that are closest to the one time point t stored in the FIFO sequence 20. In other words, δ may be from 1 to d_(ij)−1. In this case, the calculating section 140 may use the second weights w_(ij) ^([δ]) corresponding respectively to the closest d_(ij)−1 input values x_(i) ^([t−δ]).

In this case, the learning processing section 150 may update each second weight w_(ij) ^([δ]) as shown in the following expression.

$\begin{matrix} \left. w_{ij}^{\lbrack\delta\rbrack}\leftarrow{w_{ij}^{\lbrack\delta\rbrack} + {\eta\frac{1}{\sigma^{2}}{x_{j}^{\lbrack{t - \delta}\rbrack}\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)}}} \right. & {{Expression}\mspace{14mu} 25} \end{matrix}$

The calculation performed by the calculating section 140 using the closest input values and the corresponding second weights w_(ij) ^([δ]) is not limited to the calculation of the node value B_(j) ^([t]). The calculating section 140 may use the closest input values and the corresponding second weights w_(ij) ^([δ]) in the same manner to calculate the conditional probability according to Expression 4. In this way, when calculating the node value B_(j) ^([t]) and the conditional probability, the calculating section 140 can suitably reflect each of the closest input values believed to have a greater affect, and can therefore increase the learning power.

The calculation of the node value B_(j) ^([t]) performed by the calculating section 140 may include a third weight instead of the second weight. The third weight may be a value μ_(l) ^([t−δ])x_(i) ^([t−δ]) for weighting a value obtained by amplifying or attenuating, with a second coefficient μ_(l) according to the time point, the input values at each of the two or more time points between the one time point t and the time point that is a first number d_(ij) of time points before the one time point t. Specifically, the calculating section 140 may calculate the node value B_(j) ^([t]) using the d_(ij)−1 input values x_(i) ^([t−δ]) that are closest to the one time point t stored in the FIFO sequence 20. In other words, δ may be from 1 to d_(ij)−1, as shown in the expression below.

$\begin{matrix} {B_{j}^{\lbrack t\rbrack} = {b_{j} + {\sum\limits_{i = 1}^{I + H}\left( {{\sum\limits_{k \in K}^{\;}{u_{i,j,k}\alpha_{i,j,k}}} + {\sum\limits_{\delta = 1}^{d_{ij} - 1}{\mu_{i}^{\lbrack\delta\rbrack}x_{i}^{\lbrack{t - \delta}\rbrack}}} - {\sum\limits_{l \in L}^{\;}{v_{j,i,l}\gamma_{i,l}}}} \right)}}} & {{Expression}\mspace{14mu} 26} \end{matrix}$

The above describes an example in which the learning apparatus 100 of the second modification updates the node value B_(j) ^([t]) using the variance parameter σ indicating the variance in the probability distribution of the input values. However, the update of the node value B_(j) ^([t]) performed by the learning apparatus 100 is not limited to this. The learning apparatus 100 may be configured to update the node value B_(j) ^([t]) by using σ² instead of the variance parameter σ, and by using a natural gradient or the like, as shown in the expression below.

$\begin{matrix} \left. \sigma^{2}\leftarrow{\sigma^{2} + {\eta_{2}{\sum\limits_{j = 1}^{I}\left( {\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)^{2} - \sigma^{2}} \right)}}} \right. & {{Expression}\mspace{14mu} 7} \end{matrix}$

Furthermore, the learning apparatus 100 of the second modification can determine the node value B_(j) ^([t]) and the corresponding variance parameter σ by learning, and therefore may be configured to further perform a function using this variance parameter σ. For example, the learning apparatus 100 may be configured to detect peculiarity of the time-series input data of a learning target by using a learned model corresponding to the time-series input data of the investigation target.

In this case, the acquiring section 110 may be configured to acquire time-series input data of an investigation target that is a time-series of input data including a plurality of input values and different from the time-series input data used when performing the learning. The supplying section 120 may supply each input node of the model 10 with the time-series input data of the investigation target. The calculating section 140 may be configured to calculate the node value B_(j) ^([t]) of a first node corresponding to the time-series input data of the investigation target by using each propagation value propagated to the first node. Furthermore, the learning processing section 150 may update each parameter.

In this way, the learning apparatus 100 can calculate the variance parameter σ and the node value B_(j) ^([t]) of a first node corresponding to the time-series input data of the investigation target. The learning apparatus 100 may then calculate the peculiarity of the input value corresponding to the first node of the time-series input data of the investigation target by comparing the variance parameter σ_(m) corresponding to the time-series input data of the investigation target to a variance parameter σ₀ learned using the time-series input data of a learning target.

For example, if the expression shown below is established, the learning apparatus 100 may determine an input value x_(j) ^([t]) to be peculiar. The constant C₀ may be a predetermined value. The constant C₀ is 3, for example. Furthermore, the learning apparatus 100 may output, as the peculiarity of an input value, ε shown in Expression 27. In this way, the learning apparatus 100 can easily detect whether the data of the investigation target input in time series includes peculiar data that falls outside a range predicted by the learning. ε=|x _(j) ^([t]) −B _(j) ^([t]) |>C ₀σ₀  Expression 28:

Furthermore, the learning apparatus 100 of the second modification can determine the node value B_(j) ^([t]) and the corresponding variance parameter σ by learning, and therefore may be configured to generate new time-series data. The learning apparatus 100 may generate new time-series data based on a probability distribution having a variance corresponding to the variance parameter σ and having the average value of the input values at one time point as the node value B_(j) ^([t]). The learning apparatus 100 can easily generate new time-series data corresponding to the learning by using random numbers or the like.

FIG. 10 shows an example of a computer 800 in which aspects of the present invention may be wholly or partly embodied. A program that is installed in the computer 800 can cause the computer 800 to function as or perform operations associated with apparatuses of the embodiments of the present invention or one or more sections (including modules, components, elements, etc.) thereof, and/or cause the computer 800 to perform processes of the embodiments of the present invention or steps thereof. Such a program may be executed by the CPU 800-12 to cause the computer 800 to perform certain operations associated with some or all of the blocks of flowcharts and block diagrams described herein.

The computer 800, according to the present embodiment, includes a CPU 800-12, a RAM 800-14, a graphics controller 800-16, and a display device 800-18, which are mutually connected by a host controller 800-10. The computer 800 also includes input/output units such as a communication interface 800-22, a hard disk drive 800-24, a DVD-ROM drive 800-26 and an IC card drive, which are connected to the host controller 800-10 via an input/output controller 800-20. The computer also includes legacy input/output units such as a ROM 800-30 and a keyboard 800-42, which are connected to the input/output controller 800-20 through an input/output chip 800-40.

The CPU 800-12 operates according to programs stored in the ROM 800-30 and the RAM 800-14, thereby controlling each unit. The graphics controller 800-16 obtains image data generated by the CPU 800-12 on a frame buffer or the like provided in the RAM 800-14 or in itself, and causes the image data to be displayed on the display device 800-18.

The communication interface 800-22 communicates with other electronic devices via a network 800-50. The hard disk drive 800-24 stores programs and data used by the CPU 800-12 within the computer 800. The DVD-ROM drive 800-26 reads the programs or the data from the DVD-ROM 800-01, and provides the hard disk drive 800-24 with the programs or the data via the RAM 800-14. The IC card drive reads programs and data from an IC card, and/or writes programs and data into the IC card.

The ROM 800-30 stores therein a boot program or the like executed by the computer 800 at the time of activation, and/or a program depending on the hardware of the computer 800. The input/output chip 800-40 may also connect various input/output units via a parallel port, a serial port, a keyboard port, a mouse port, and the like to the input/output controller 800-20.

A program is provided by computer readable media, such as the DVD-ROM 800-01 or the IC card. The program is read from the computer readable media, installed into the hard disk drive 800-24, RAM 800-14, or ROM 800-30, which are also examples of computer readable media, and executed by the CPU 800-12. The information processing described in these programs is read into the computer 800, resulting in cooperation between a program and the above-mentioned various types of hardware resources. An apparatus or method may be constituted by realizing the operation or processing of information in accordance with the usage of the computer 800.

For example, when communication is performed between the computer 800 and an external device, the CPU 800-12 may execute a communication program loaded onto the RAM 800-14 to instruct communication processing to the communication interface 800-22, based on the processing described in the communication program. The communication interface 800-22, under control of the CPU 800-12, reads transmission data stored on a transmission buffering region provided in a recording medium such as the RAM 800-14, the hard disk drive 800-24, the DVD-ROM 800-01, or the IC card, and transmits the read transmission data to network 800-50 or writes reception data received from network 800-50 to a reception buffering region or the like provided on the recording medium.

In addition, the CPU 800-12 may cause all or a necessary portion of a file or a database to be read into the RAM 800-14, the file or the database having been stored in an external recording medium such as the hard disk drive 800-24, the DVD-ROM drive 800-26 (DVD-ROM 800-01), the IC card, etc., and perform various types of processing on the data on the RAM 800-14. The CPU 800-12 may then write back the processed data to the external recording medium.

Various types of information, such as various types of programs, data, tables, and databases, may be stored in the recording medium to undergo information processing. The CPU 800-12 may perform various types of processing on the data read from the RAM 800-14, which includes various types of operations, processing of information, condition judging, conditional branch, unconditional branch, search/replace of information, etc., as described throughout this disclosure and designated by an instruction sequence of programs, and writes the result back to the RAM 800-14. In addition, the CPU 800-12 may search for information in a file, a database, etc., in the recording medium. For example, when a plurality of entries, each having an attribute value of a first attribute is associated with an attribute value of a second attribute, are stored in the recording medium, the CPU 800-12 may search for an entry matching the condition whose attribute value of the first attribute is designated, from among the plurality of entries, and reads the attribute value of the second attribute stored in the entry, thereby obtaining the attribute value of the second attribute associated with the first attribute satisfying the predetermined condition.

The above-explained program or software modules may be stored in the computer readable media on or near the computer 800. In addition, a recording medium, such as a hard disk or a RAM provided in a server system connected to a dedicated communication network or the Internet, can be used as the computer readable media, thereby providing the program to the computer 800 via the network.

The present invention may be a system, a method, and/or a computer program product. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.

The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.

Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.

Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like, and conventional procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to individualize the electronic circuitry, in order to perform aspects of the present invention.

Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.

These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.

The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.

The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.

While the embodiments of the present invention have been described, the technical scope of the invention is not limited to the above described embodiments. It is apparent to persons skilled in the art that various alterations and improvements can be added to the above-described embodiments. It is also apparent from the scope of the claims that the embodiments added with such alterations or improvements can be included in the technical scope of the invention.

The operations, procedures, steps, and stages of each process performed by an apparatus, system, program, and method shown in the claims, embodiments, or diagrams can be performed in any order as long as the order is not indicated by “prior to,” “before,” or the like and as long as the output from a previous process is not used in a later process. Even if the process flow is described using phrases such as “first” or “next” in the claims, embodiments, or diagrams, it does not necessarily mean that the process must be performed in this order.

As made clear from the above, with the embodiments of the present invention it is possible to learn a model corresponding to real number time-series input data and to increase the expressive power and learning power of a learning apparatus. 

What is claimed is:
 1. A learning method for learning a model corresponding to time-series input data, comprising: acquiring the time-series input data, the time-series input data being a time series of input data including a plurality of input values having non-binary values; propagating, to a plurality of hidden nodes in the model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point; generating update parameters of the plurality of hidden nodes based on the plurality of input values, at least one of the update parameters having a non-binary value; calculating a node value of a first node among the plurality of hidden nodes using each propagated value propagated to the first node, the node value including an average value of possible values of the first node; and updating, using a processor, a weight parameter used to calculate each propagation value propagated to the first node using a corresponding input value, a variance in a probability distribution of input values, and a calculated error of the node value at the one time point.
 2. The learning method according to claim 1, further comprising updating, for the first node, a variance parameter indicating the variance in the probability distribution of an input value, by using the corresponding input value and the calculated error of the node value at the one time point.
 3. The learning method according to claim 2, wherein the variance parameter is a parameter common to the plurality of input values.
 4. The learning method according to claim 3, wherein updating the weight parameter includes making an update amount of the weight parameter smaller when the variance parameter is larger.
 5. The learning method according to claim 3, wherein updating the variance parameter includes updating the variance parameter based on a sum of mathematical squares of the calculated error for each of the plurality of hidden nodes.
 6. The learning method according to claim 1, wherein a portion of the plurality of input values includes multi-value input values.
 7. The learning method according to claim 1, wherein the plurality of input values includes multi-value input values and binary input values.
 8. The learning method according to claim 7, further comprising, for a second node associated with an input value that is a binary value among the plurality of input values, calculating a conditional probability of the second node being one of two values of the binary value, based on each propagation value propagated to the second node and an input value at the one time point corresponding to the second node.
 9. The learning method according to claim 7, further comprising: for a second node that is associated with an input value that is a binary value among the plurality of input values, calculating an expected value of a node value of the second node, based on each propagation value propagated to the second node and an input value at the one time point corresponding to the second node; and updating a weight parameter used to calculate each propagation value propagated to the second node, by using the input value corresponding to the one time point and an error of the expected value of the node value of the second node.
 10. The learning method according to claim 1, further comprising: updating a first update parameter of the update parameters based on an input value at a time point that is a first number of time points before the one time point; and calculating the node value of the first node based on a value obtained by weighting the first update parameter with a first weight included in the weight parameter.
 11. The learning method according to claim 10, wherein updating the first update parameter includes updating the first update parameter by multiplying a first coefficient by a sum of original first update parameters and the input value at the time point that is the first number of time points before the one time point.
 12. The learning method according to claim 10, wherein the weight parameter includes a second weight corresponding to input values at each of two or more time points between the one time point and the time point that is the first number of time points before the one time point.
 13. The learning method according to claim 10, wherein the weight parameter includes a third weight for weighting a value obtained by amplifying or attenuating, with a second coefficient corresponding to time points, the input values at each of two or more time points between the one time point and the time point that is the first number of time points before the one time point.
 14. A learning method using a learned model corresponding to time-series input data of a learning target, comprising: acquiring time-series input data of an investigation target, the time-series input data being a time series of input data including a plurality of input values having non-binary values; propagating, to a plurality of hidden nodes in the learned model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point; generating update parameters of the plurality of hidden nodes based on the plurality of input values, at least one of the update parameters having a non-binary value; calculating a node value of a first node among the plurality of hidden nodes using each propagated value propagated to the first node, the node value being including an average value of possible values of the first node; and detecting, using a processor, peculiarity of an input value corresponding to the first node and input at the one time point, using the node value of the first node and a variance parameter indicating a variance in a probability distribution of input values learned using the time-series input data of the learning target.
 15. The learning method according to claim 14, wherein at least a portion of the plurality of input values includes multi-value input values.
 16. The learning method of claim 14, further comprising: updating a first update parameter of the update parameters based on an input value at a time point that is a first number of time points before the one time point; and calculating the node value of the first node based on a value obtained by weighting the first update parameter with a first weight included in the weight parameter.
 17. A non-transitory computer readable storage medium having instructions embodied therewith, the instructions executable by a processor or programmable circuitry to cause the processor or programmable circuitry to perform operations comprising: acquiring the time-series input data, the time-series input data being a time series of input data including a plurality of input values having non-binary values; propagating, to a plurality of hidden nodes in the model, each of a plurality of propagation values obtained by weighting each input value at a plurality of time points before one time point according to passage of time points, in association with the plurality of input values at the one time point; generating update parameters of the plurality of hidden nodes based on the plurality of input values, at least one of the update parameters having a non-binary value; calculating a node value of a first node among the plurality of hidden nodes using each propagated value propagated to the first node, the node value being including an average value of possible values of the first node; and updating, using a processor, a weight parameter used to calculate each propagation value propagated to the first node, using a corresponding input value, a variance in a probability distribution of input values, and a calculated error of the node value at the one time point.
 18. A non-transitory computer readable storage medium having instructions embodied therewith, the instructions executable by a processor or programmable circuitry to cause the processor or programmable circuitry to perform the learning method of claim
 14. 19. The method of claim 1, wherein the probability distribution of the input values is expressed as: ${p\left( x_{j}^{\lbrack t\rbrack} \right)} = {\frac{1}{\sqrt{2\pi}\sigma}\exp\left\{ {- \frac{\left( {x_{j}^{\lbrack t\rbrack} - B_{j}^{\lbrack t\rbrack}} \right)^{2}}{2\sigma^{2}}} \right\}}$ where x_(j) ^([t]) is an input value to a j^(th) node at a time t, σ is the variance, and B_(j) ^([t]) is the average value of possible values of the j^(th) node at time t. 